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ABSTRACT 


We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of 
the planetary nebula NGC 7009. The primary goal was to obtain high spatial 
resolution of the intrinsic line ratio [O Til] 4364/5008 and theicby evaluate tin. 
electron temperature (T e ) and the fractional mean-square T e variation (t\) across 
the nebula. The WFPC2 T e map is rather uniform; almost all values are between 
9000 1 1,000 K, with the higher T e s closely coinciding with the inner He ++ -zone. 

The results indicate very small values - <0.01 - for i l A throughout. Our STIS 
data allow an even more direct determination of T e and t A , albeit for a much 
smaller area than with WFPC2. VVe present results from binning the data along 
the slit into tiles that are 0.5" square (matching the slit width). The average 
[() in] temperature using 45 tiles (excluding the central star and STIS fiducial 
bars) is 10,139 K; t\ is 0.0035. The remarkably low values measured for t\ 
are a significant new find. However, these pertain to only 2-dimensions. The 
observations here do not address T e fluctuation along the line of sight. Because 
each element of area treated in the plane of the sky represents a column which 
has already created a spatially averaged temperature along the line of sight, it is 
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likely that the value for the fractional mean-square T e variation in 3-dimensions 
(ty) is substantially higher than i\. Measurements of the [fractional] mean- 
square T e variation along the line of sight {t 2 ) are needed to reliably measure ty. 
Even with our determination of very low values for t\, we cannot rule out the 
possibility that temperature fluctuations may explain for NGC 7009 the large 
discrepancy between heavy element abundances inferred from emission lines that 
are eollisionallv excited compared with those that are due to recombination lines. 

Subject headings: ISM: planetary nebulae — ISM: abundances ISM: atoms - 
ISM: individual (NGC 7009) 


1. INTRODUCTION 

Most observational tests of the chemical evolution of the universe rest on emission line 
objects; these define the endpoints of stellar evolution and probe the current state of the 
interstellar medium. Gaseous nebulae are laboratories for understanding physical processes 
in all emission-line sources, and probes for stellar, galactic, and primoidial nucleosynthesis. 

There is a fundamental issue that continues to be problematic the discrepancy be- 
tween heavy element abundances inferred from emission lines that are eollisionallv excited 
compared with those due to recombination/cascading, the so-called iecombination lines . 
Studies of planetary nebulae (PNs) contrasting recombination and collisional abundances 
(Liu et al. 1995, K witter k Henry 1998) often find differences exceeding a factor of two. In 
an extensive study of NGC 6153, Liu et al. (2000) found that C++/H+, N ++ /H + , 0++/H + , 
and Ne ++ /H + ratios derived from optical recombination lines are all a factor of ~10 higher 
than the corresponding values deduced from collisionally-excited lines. Abundances deter- 
mined from these two methods disagree by a factor larger than the spread of abundances 
used to determine such fundamental quantities as Galactic abundance gradients (e.g., Shaver 
et al. 1983; Simpson et al. 1995; Henry k Worthey 1999). 

Most of the efforts to explain the abundance puzzle between collisional and recombi- 
nation values have attempted to do so by examining electron temperature (7),) variations 
in the plasma. This is often done, using the formalism of Peimbert (1967), in terms of the 
[fractional] mean-square variation (t 2 ) of T e . The inferred metallicity obtained by using the 
usual (optical/UV) forbidden lines is very sensitive to T e (exponential) and t 2 . On the other 
hand, recombination lines are rather insensitive to T e and t 2 . Agreement, close to the higher 
recombination value can be forced in the derived abundance by attributing the difference to 
(solving for) I 2 . Lot us consider the case of the PN NGC 7009, also known as the “Sat urn 
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Nebula”; its image produced by Balick et al. (1998) from WFPC2 data has graced the cover 
of coffee-table books (Petersen k Brandt 1998). NGC 7009 has stood near the centre of 
the abundances controversy. It has long been known that its recombination lines are very 
strong, particularly those of 0 IT. Liu et al. (1995) found that the recombination C, N, and 
() abundances are a factor of —5 larger than the corresponding collisional abundances. This 
was recently found to be the case also for neon (Luo, Liu, k Barlow 2001). One possible 
reconciliation of these two drastically different abundances, which Liu et al. (1995) discuss, 
is to invoke large I 2 . For oxygen in NGC 7009, by invoking t 2 ~ 0.1, agreement can be forced 
in the abundance derived from optical lines close to the higher recombination value a value 
more than 2.5 times larger than the solar O/H of 7.41 x HO 4 (Grevesse k Sauval 1998). Such 
a large t 2 is not at all predicted by current theory/models (e.g., Kingdon k Ferland 1998). 

The current unsettled situation has led to efforts to broaden the study to include othei 
variables besides T r to analyse the effects upon abundance determinations. One promising 
avenue is to examine this considering density variations, abundance variations, and 7 ( , wui- 
ations in combination. Liu et al. (2000) took this approach in their investigation of the PN 
NGC 0153 with a two-phase empirical model. Recall this object, has even more disparate 
abundance differences (factor of —10) than does NGC 7009. Pequignot et al. (2002) contin- 
ued the study using photoionization models including two components with different heavy 
element abundances. Basically, one phase is small, dense, ionized clumps that aie highly 
enhanced in heavy (dements and deficient in hydrogen. This phase would have T e — 10' 1 K 
and contribute almost all the heavy element recombination line emission. These clumps aie 
embedded in the second phase which has the more usual nebular properties/composition, 
with T t ~I0 4 K, and which is responsible for almost all the emission in the optical/UV 
eollisionally excited lines. The latter phase comprises the great bulk of the nebular mass. 
There is ongoing work by this same group to try to explain other PNs of this kind, including 
the Galactic bulge PN M 1-42, which has a factor of —20 abundance dichotomy (Liu et al. 
2001) and Hf 2-2. which has the most extreme abundance difference to date, a factor of ~80 
(Liu 2002). 

In this paper, we focus solely on T e variations with the purpose to test whether the 
observational data support the possibility that t 2 in NGC 7009 may be as large as -0.1. 
In section 2, we present the HST observations and data reduction procedures. Section 3 
describes our techniques for obtaining emission line fluxes from the YVFPC2 data and how 
we use cospatial STIS data to help this process. In section 4, there is a discussion and 
analyses of extinction. YVe determine the electron temperature distributions in section 5. In 
section (i, we analyse the T e distributions in terms of average temperatures and fractional 
mean-square temperature variations in the plane of the sky. Section 7 provides a discussion 
and conclusions. 
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2. HST OBSERVATIONS AND DATA REDUCTION 

The observations of NGC 7009 described here were taken as part of our HST Cycle 8 
program GO-8114. For both the VVFPC2 and STIS visits, we set our positions with respect 
to the position of the central star (CS): a, 8 = 21 h 04 m 10>!81, -11°21'47"91) (equinox .12000). 


2.1. WFPC2 Observations 

VVe observed with VVFPC2 on 2000 April 6-7 (UT) as the first visit of the GO-8114 
program. Observations (with corresponding line of interest) of NGC 7009 were made in line 
filters F437N ([() 111] 4364 A), F487N (H/l 4863 A), F502N ([O III] 5008 A), and F656N 
(Ha 6565 A), plus continuum filter F547M. All wavelengths in this paper are vacuum rest 
wavelengths. The images were placed in the WF3 chip (with the CS at WF3-FIX). The pixel 
size is 0.0996". In order to clean cosmic rays, two images were taken in each filter except, for 
three images in F437N. Total integration times were respectively: 2200, 520, 320, 400, and 
240 s. 

The data sets we processed (also for STIS) were those obtained after sufficient time 
elapsed from the observation dates in order that “best reference files” would be stable/ finalized 
VVe requested On-the-Fly Calibration for science files and Best Reference Files. VVe then co- 
added and cosmic ray cleaned images using standard packages in IRAF. 10 Flux calibration 
for each filter was done by basically following the procedures used for the VVFPC2 Exposure 
Time Calculator. These account for the integration time, analog-to-digital converter gain of 
7 (used for F437N and F487N) or 14 electrons/DN (used for F502N, F547M, and F656N), 
and the total system throughput. 

The calculations for the contributions to the narrow-band filter observations of the 
continuum and other lines will be discussed later (§3) after the STIS observations section. 


2.2. STIS Observations 

On 2000 August 9 (UT), STIS long-slit spectra of NGC 7009 were taken. These data 
comprise 3 orbits and were taken at position angle (PA) = 86.72°. Our original plan was to 
align the slit with the major axis of NGC 7009, which is at PA = 78.5°. However, we were 
advised by STScI that there were no guide star pairs for our planned visits. VVe settled on 


1(1 [R AF is distributed by NOAO, which is operated by AURA, under cooperative agreement with NSK. 
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the revised PA, a change of only 8.22°, as this was the closest to the major axis that allowed 
for a pair of guide stars. We used a slit width of 0.5" and a slit length of 52" for these CCD 
observations. 

Because the central star (CS) in NGC 7009 is reasonably bright (V = 12.78, B = 12.66), 
it would have caused problems (bleeding, saturation, etc.), especially on long exposures, if 
it were in our science slit. On the other hand, the CS provides an excellent fiducial point in 
the nebula, where we may assess the positional registration of emission lines in oui vaiious 
exposures with STIS. Thus, we observed with the CS off the upper edge of our 52" x 0.5" 
slit by 0.1". In other words, our slit is defined with the CS moved 0.35" (0.25 + 0.1) in the 
dispersion direction (~NNW). Furthermore the slit was not centred on the CS in the spatial 
direction. Because the structure on the west-side of the nebula was of more interest to us, we 
displaced the slit centre 3.5" from the CS in the spatial direction (~W along PA = 86.72°). 
Spectra were taken with gratings G430M (with wavelength settings: 3680. 4194, 4451, 4706, 
4961) and G750M (settings: 5734, 6581). Each exposure was done in accumulation mode 
and at least two spectra were taken at each setting in order to cosmic-ray (CR) clean. After 
retrieving the data sets, we then co-added and cosmic ray cleaned images using standard 
packages in 1RAF. Calibrations to produce 2-dimensional (2D) rectified images were then 
carried out. From these, we singled out specific emission lines for further investigation. 

Data for the [() Til] 4364 A line as well as Hy 4342 A were contained in the G430M/4451 
grating setting (150s); for the [O III] 5008, 4960 A lines and H/f 4863 A lines, G430M/4961 
was used (80s); for Ha 6565 A and the [N II] 6550, 6585 A lines, G750M/6581 was used 
(90s). For G430M, the dispersion is 0.28 A per pixel for a point source and the plate scale 
0.05 arcsec/ pixel (Leitherer et, al. 2001, Chapter 13). For a uniformly filled slit with width 
0.5", a degradation in resolving power by a factor of 10 is expected to a spectral resolution 
of 2.8 A. For G750M, the dispersion is 0.56 A per pixel for a point source and the plate scale 
0.05 arcsec/pixe). For a uniformly filled slit with width 0.5", the spectral resolution would be 
5.6 A. For the analysis presented here, we were interested mainly in the distribution of line 
llux along the slit spatial direction. This was accomplished with the IRAF routine blkavy in 
conjunction with specialized software tools that we developed omselvcs. 

Even after applying the standard CR rejection there still remain many bad pixels due 
to CRs and/or hot pixels. There is considerable danger that including these can corrupt the 
llux values we seek. The program developed to eliminate these remaining bad pixels is called 
PIXHUNTER, which is described briefly in Appendix A. Once the columns containing the 
line have been cleaned for deviant pixels, we are ready to subtract an equivalent spectral 
range of continuum. We do this by using IRAF functions, including blkavi), to operate upon 
the appropriate sections of cleaned continuum. The 1-D distribution of line flux vs. spatial 



- 6 - 


coordinate for the various emission lines of interest is vvliat we need for our subsequent 
physical analysis. 

We note that there is excellent agreement with a cross check of the 1-D results of flux 
vs. spatial direction by comparing with 1-D results of flux vs. wavelength for a corresponding 
spatial sample. The latter were measured with the splut package. With this, the underlying 
continuum is fitted and the integrated line flux determined with the e-option (area under 
the line profile), which was preferable to fitting with a Gaussian profile. Because of the 
spectral impurity introduced by the relatively wide slit used, the line profiles have flatter 
tops and less extended bases (i.e., they are more “trapezoidal”) than the Gaussian fits. It is 
also apparent that the Gaussian fit is overestimating the line flux. 

Both the [O ITT] 5008 and 4960 A lines were observed simultaneously with the G430M/4961 
grating setting. We desired to accurately measure the ratio of these line fluxes to compaie 
with the most recent calculation for the theoretical intensity ratio 2.984 (Storey & Zeippen 
2000). Because both transitions arise from the same upper level, the intrinsic flux ratio 
depends only on the transition probabilities (A-values) and wavelengths. What we found 
was a surprising variation in the F(5008)/F(4960) ratio with position along the slit. 1 his 
amounts to a variation in the ratio of roughly 3.0±0.1. Furthermore, it appears more-or-less 
periodic with an cycle. 

According to Ted Gull (private communication), thus is an instrumental effect and is a 
ratio of two fringe patterns. The source of the problem is a thin blocker filter that had to 
be matched with each grating and the best (and oidy) location that it could be placed was 
above the grating in a stable mounting. There is some discussion of fringing at the S I IS 
web site (http://www.stsci.edu/instruments/stis/performance/anomalies). More specific in- 
formation appears in chapter 7 of the STIS Instrument Handbook (Leitherer et. al. 2001) 
under a section called “Fringing due to the Order Sorter filters”. See 

(www.stsci.edu/instruments/stis/documents/liandbooks/cyclel l/c07_performance2.html#326437). 

This distinct fringing pattern also appears in similar STIS data of the Orion Nebula obtained 
under program GO-7514 (PI HR). 

To attempt, to do anything about fringing would probably require a dedicated HST/STIS 
calibration program. We do not treat this further here. There is no apparent evidence 
of fringing in the F(5008)/F(4364) ratio. As will be seen, the results for t he temperature 
analyses are in good accord between the two datasets using STIS and WFPG2, which provides 
some support that any fringing in the STIS data is not affecting our conclusions. Because 
of the fringe pattern, we measured the fluxes F(5008) anti F(4960) over an integer number 
of fringing cycles and then calculated the F(5008)/F(4960) ratio. The overall value of 3.008 
is in good agreement with the theoretical prediction, given that extinct ion will raise slightly 



- 7 

the observed ratio above the intrinsic one (see §4). 


3. EXTRACTION of EMISSION LINE FLUXES from WFPC2 DATA 

One of the most important operations involved partitioning the contributions to the 
F437N emission. We consider these as due to [O III] 4364, H'y, and continuum emission. 
(An examination of our STIS spectra shows no other emission lines will contribute significant 
flux.) It is necessary to take into account the total system throughput for the F437N filter 
(Biretta et al. 2001). A finer scale digital version of the total system throughput for this 
or any WFPC2 filter can also be produced using the task calcband. We used a digital 
form to obtain the throughputs for both the [O III] 4364 and Hy lines accounting for the 
Doppler correction for the actual velocity of NGC 7009 at the time of obseivation. 1 his was 
—20.71 km s _1 after accounting for the heliocentric velocity of 7009, which is —46.6 km s 
(Schneider 1983). While still a net blue shift (-0.3 A), this does result in a slightly higher 
H'y leakage due to near maximum red shift of Earth’s motion at the time of obseivation. 
The relative transmission at the Doppler corrected wavelengths for Hy to 4364 was 0.0253. 
Because the transmission of F437N for the Hy line is much less than for the 4364 line, this 
is termed Hy “leakage” . 

In order to assess the contribution of 4364 and Hy leakage to our F437N data, we did a 
careful spatial registration/comparison of the highest signal-to-noise ( S / N) portions of our 
STIS long-slit, data with corresponding WFPC2 filters. Two areas were chosen that we call 
El and WW that were regions of higher S/N and avoided the CS as well as the fiducial bars. 
They span from 1.05" to 7.60" east and from 0.75" to 8.10" west of the CS, respectively, 
along our 0.5"-wide STIS slit. The registration is very close with the El area from STIS 
3.275 and from WFPC2 3.2737 arcsec 2 . For WW, areas are respectively, 3.675 and 3.6705 
a reset: 2 . 

Initially, we also tried to determine the continuum contribution to F437N from the 
spatial registration of our STIS data with the corresponding WFPC2 filters. All our efforts 
resulted in continuum fluxes that were unreliable and generally too large for what could 
possibly fit for the WFPC2 observations. While it is not clear why the flux levels underlying 
the emission lines (not just 4364 but also Ha, H/J, and 5008) were invariably too high, the 
integration times with STIS were insufficient to provide adequate signal-to-noise to measure 
these continuum fluxes. 

Because t he dominant contribution to the nebular continuum for NGC 7009 is recombi- 
nation processes, we estimate the continuum contribution to each of the WFPC2 filters from 
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recombination theory. These were provided independently by two of us (JPH and G.JF ) from 
theoretical continuum models, with similar results. For these, we assumed T e = 10,000 K, 
N(He)/N(H) = 0.11, and ionic fractions: N(He + )/N(He) = 0.88 and N(He ++ )/N(He) = 
0.12, which are essentially the NGC 7009 values in Kingsburgh & Barlow (1994). The scaled 
values adopted for the continuum emission relative to the H/f line emission (all intrinsic) are 
6.91E-4, 6.53E-4, 6.45E-4, and 5.65E-4 A -1 for F437N, F487N, F502N, and F656N, respec- 
tively. To obtain the total continuum contribution in the respective filters, one multiplies 
by the appropriate “rectwidth” (RECTW) 31.994, 33.8, 35.502, and 29.89 A obtained using 
the task bandpar . 1 1 

This provides directly the predicted continuum correction for F487N of 0.02207 l( H/i), 
where I(H/J) is the intrinsic flux. Because both line and continuum are equally affected by 
extinction, we may also work with the observed fluxes without regard to extinction correction 
here. To predict the continuum correction for F656N and F437N, we scale the H/f result to 
the Ha and H" lines using recombination theory (Storey & Hummer 1 99;)) assuming 7„ = 
10,000 K and A r e = 10,000 cm -3 , Case B: I(Ha)/I(H/?) = 2.85 and l(H/i)/l(H'>) = 2.13. 
This results in predicted continuum corrections for F656N of 0.005920 I(Ho) and for F437N 
of 0.04709 I(H7). The continuum relative to 4364 is obtained by multiplying the ratio of 
STIS observed fluxes of Ht/ 4364 by the continuum/Hy value of 0.04709. The differential 
extinction from 4304 to 4342 A is negligible for 7009. 

The continuum correction for F502N relative to [O III] 5008 is obtained by multiplying 
the ratio of STIS observed fluxes of H/J/5008 by the continuum/H/J value of 0.02294. Because 
the 5008 line is the strongest of the set, the continuum correction is extremely small; hence, 
it is safe to neglect the very small differential extinction from 4863 to 5008 A. 

While the above discussion is appropriate for correcting for the continuum in the specific 
WFPC2 filters, if comparison is made between individual filters, then an extinction correction 
would need to be considered (see §4). 

Our original observing plan had been to use the WFPC2/F547M data to help determine 
the continuum emission. This does not appear to be a reliable method as found from a 
detailed comparison of our STIS and WFPC2 data. For instance, let us consider area El. 
The observed F547M flux is 1.128xl0~ 12 erg cm -2 s~ l . Using the continuum predictions of 
the theoretical model mentioned earlier for the mean wavelength 5468 A of this filter, the 
continuum emission at 5468 A relative to the H (i line emission is 6.25E-4 A -1 . To obtain the 


"Note that RECTW is not the “effective width of the bandpass” tabulated in Table G.l and also shown 
as dashed rectangles in the Passband Plots of Filter Throughput (Appendix 1) of the WFPC2 Instrument 
Handbook (BireUaet al. 2001). 
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total continuum contribution in this filter, we multiply by the R.ECTYV — (5.18. 58 A and then 
multiply by the observed STIS Hf) flux of 1.507xl0~ 12 erg cnr 2 s" 1 to arrive at G.014xl0- 13 
erg cm -2 s _1 . This is only 53% of the observed F547M flux. In the above assessment, we 
have neglected the small differential extinction. Presumably the other ~47% is due to line 
emission. The zero- intensity width of the F547M filter extends from ~5000-6000 A. 

In addition to accounting for the continuum in each YVFPC‘2 filter and the Hy leakage 
in F437N, it is also necessary to account for the contribution of [N IT] 6585.23 and 6549.86 A 
emission into F656N. This is done by measuring the fluxes of these lines, as well as that 
of Ha, in the areas El and YVYV. The method follows from that used for the Hy leakage 
in F437N. The Doppler corrected wavelengths for the three lines are determined and the 
relative throughputs applied for the F656N filter. These are 0.03645 and 0.4552 for the 6585 
and 6550 A lines respectively relative to Her. Note that although the 6585 flux will always 
be larger than that of 6550, the latter contributes more to the observed F656N flux. 

Table 1 shows the details of the comparison with the STIS overlay for El and YVYV 
for each of the YVFPC2 narrow-band filters. Column 2 has the ratios of the observed line 
flux with STIS to the corresponding narrow band filter flux with YVFPC2. Column 3 shows 
the flux ratio Contm. /emission line, derived from recombination theory and the YVFPC2 
filter profile. Columns 4 and 5 provide the observed line flux ratio derived from STIS data as 
modified passing through the pertinent YVFPC2 filter profile. Column (i then shows the Total 
(of the previous columns 2 up to 5) for the STIS flux (including the theoretical continuum) 
divided by the F437N flux. The results for both areas are in close agreement; thus we use the 
combined area El 4- YVYV, which we call SUM, for our actual calibiation to obtain obseived 
line fluxes from the YVFPC2 data. In the case of F437N and F502N, the total or the STIS 
emission plus inferred continuum exceeds the YVFPC2 flux by a factor of 1.089 and 1.043 
respectively, while for F487N and F656N, the total STIS emission plus continuum is less 
than the YVFPC2 flux by factors of 1.032 and 1.037 respectively. From the SUM entry line 
for F437N in Table 1, we see that the relative strengths of the components comprising F437N 
are roughly in the proportions 0.71 : 0.19 : 0.10 for the 4364 line : continuum : Hy leakage. 

YVe use the results for the SUM area to define our best relationships for determining 
the necessary emission line fluxes from the YVFPC2 data. Because the f 43 7N collections 
for the continuum and Hy leakage are tied to a measurement of a Balmer line flux, we use 
H/j from the F487N image to measure its flux. Thus we arrive at equations (1) - (4). 


F(Ha) = 0.9531 F656N , 
F(H0) = 0.9482 F487N , 


( 1 ) 

( 2 ) 



F(5008) = 1.041 F502N , 

F(4364) = 1.089 F437N - 0.03218 10 _tU24 c(W) F487N. 


( 3 ) 

( 4 ) 


Equation (4) requires further description. Both the H7 leakage and continuum collections 
relative to H7 are known for the specific STIS-calibrated areas above. However, when we 
deal with other parts of the nebula observed with YVFPC2 but not by the STIS slit, there 
is no direct way of knowing the observed H7 flux. Therefore, we make these corrections in 
terms of the observed H fj flux, which may be assessed via the F487N dat a. Again, we use 
the theoretical flux ratio H/J/Hy = 2.13 but must correct for differential extinction between 
the two lines. In the presence of extinction, the observed H7 line will be weaker relative 
to the observed H/l line (than the theoretical ratio). Therefore, the contribution to h 437N 
will be less and thus the correction term should be less. The extinction is measured by the 
parameter c(H/7), which is the logarithmic extinction correction at H/i, discussed in the next 
section. The constant 0.03218 in equation (4) was determined by requiring perfect agreement 
amongst all the pertinent observed fluxes for area SUM. 


4. EXTINCTION and REDDENING CORRECTION 

Before deriving the T e distribution from either the STIS or YVFPC2 data, we first correct 
for extinction. This is calculated by comparing the observed F(H«)/F(H/1) ratio with the 
theoretical ratio l(Hu!)/I(H/l). Again, we use a value of 2.85 assuming T t , = 10,000 K and 
N e = 10,000 cm -3 , Case B (Storey k Hummer 1995). The extinction correction is done in 
terms of c( H/7), given by the relationship 

log[F{\)/F(Hp)] = loy\I {X) / 1 {H P)\ - /(A) c(/7/i), (5) 

where /(A) is the extinction curve. For the 4364, 4863, 5008, and 6565 A lines, the respective 
values for /(A) are 0.124, 0, -0.034, and -0.323 (Seaton 1979). This leads to, 

c(Hp) = 3.096 lo<j\F(Ha)/ F(H P)] - 1.408. (6) 

The correction for extinction/reddening from observed to intrinsic flux for the 4364 and 
5008 lines is then given by, 

7(4364) = F(4364) 10 1 ’ 124 c{H0) ; 7(5008) = F(5008) 10 O9fiG c{u0) . (7) 

Note that the first of these was used to derive equation (4). 

For the STIS data, we binned the pixels along the slit into tiles that are 0.5" square 
(matching the slit width). The CS and STIS fiducial bars are excluded. 1 his produced 
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r( H/j) results for 45 tiles, which are the same set used later for the T e analysis. The statistics 
without regard to any weighting for brightness are: the average c(H/J) is 0.1032; the median 
is 0.1038; and the values range from 0.0239 to 0.1592. The median and average are small 
and close to previous spectroscopic values, e.g., 0.12 (Barker 1983 and references therein). 

For the YVFPC2 data, we work with the individual pixels. There were a few negative 
r(H/i) values, which are replaced with a value of zero. There also were spuriously high c(H 0) 
values inferred in the FLIERS (Fast Low Ionization Emission Regions), particularly on the 
west side. The reason for this is that the [N II] 6550, 6585 A lines are very much stronger in 
the FLIERS. This can be seen from the paper by Balick et al. (1998) in their colour image 
(fig. 1) and emission-line profiles (fig. 3). The leakage by the 6550, 6585 A lines into F656N 
enormously exceeds the ^calibration” shown in Table 1. The application of equation (1) 
overestimates F(H«) in the FLIERS, resulting in a spuriously high F(Ha)/F(H/i) and thus 
(•(H/i), reaching as high as 0.456. With the guidance of the STIS slit overlay, we chose to cap 
the r(H/7) value in the vicinity of the FLIERS at a value of 0.14. This value is not crucial 
and could have been set between 0 - 0.2 because extinction is low for NGC 7009. As will be 
seen in the next section, it matters little for the subsequent temperature analysis. 


5. ELECTRON TEMPERATURE DETERMINATION 

The electron temperature T e is derived from the intrinsic ratio 1(5008)/I(4364) using 
the following relation, 

T e = 32,96G/[Jn(/(5008)//(4364)) - 1.70 1 ] . (8) 

Effective collision strengths are from Burke, Lennon, k Seaton (1989) for T e = 10 4 K. 
Transition probabilities (A-values) are from Froese Fischer k Saha (1985). Note that this 
holds in the low-A r f limit, which should be valid for NGC 7009 where N e values are less than 
the critical densities (Ncnt) for these lines. The lowest A r crjf ~6.4xlO J cm 3 for the 5008 
line, which is well above A r e values determined (e.g., Hyung k Aller 1995a, b). 

For the STIS data, we continue with the analysis using the 45 tiles along the slit de- 
scribed above. Equations (6)-(8) are applied to the four emission lines to derive T e . Figure 1 
shows the distribution of T e vs. tile number, which increase from ~E to ~W along the slit. 
The dashed line represents a linear interpolation across the tiles that were deemed to have 
unreliable measurements because of proximity to the CS (25-27) and the east fiducial bar 
(9 11). There appears to be some symmetry with respect to the CS. From a more or less 
flat T r ~ 10,800 K over the centre tiles, there is a decrease to ~ 9500 K and then a rise again 
to roughly 10,800 K in bot h the east, and west with increasing distance from the star. The 
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situation is less dear on the east side, partly because of the lack of data at the east fiducial 
bar. In the next section, we will evaluate this distribution in terms of T ( . variations. 

A plot of the observed relative [O III] 5008 A flux is also shown in Figure 1. This curve 
is displayed (unsmoothed) at the pixel level providing 0.05” spatial resolution along the slit. 
There are some noteworthy oscillations between ~8--13" west of the CS. There appears to 
be 4 cycles with a period between ~1 1.5" per cycle. These reach a factor of ~2 brightness 

contrast between adjacent local maxima and minima. The structure is very repeatable when 
compared to the corresponding run of observed Ha and H/f llux. Because t hese oscillations 
do not appear to be correlated with the T e distribution and are not present in the flux ratios 
of 5008 to the Balmer lines, they must be due to a column density (or j A J jdl) modulation. 
Such a pattern might result from episodic mass loss. We find very similar spatial behavior, 
including this tantalizing oscillatory pattern, from our WFPC2 F502N, F187N, and F65GN 
data when we overlay the STIS slit. We note that this result from our 1-F) cut is insufficient 
to establish a case for episodic mass loss. 

Using the WFPC2 data, the path to T e is more uncertain than is the case with our 
STIS measurements. However, the WFPC2 data cover a much larger area of the nebula. 
Equations (1) (4) and (6)-(8) are used to derive T e . In Figure 2, we present the 7’,. map 
at a resolution of 0.1". The outer region with white and black speckles is omitted from the 
analysis. The S/N in F437N is generally too low out there to obtain reliable f) information. 
Nevertheless, with a hint of the ansae, this area helps visually identify the major axis. We 
note that by binning the data in the outlying areas, the S/N may be improved. For the 
purposes of this paper, we restrict the analyses to the interior main body of the nebula with 
the highest S/N, where there is sufficient area to address temperature fluctuations in this 
nebula. 

The high T c white region immediately surrounding the CS is the result of diffracted or 
scattered stellar continuum. Equation (4), which corrects for the continuum in the F437N 
filter, does not account for stellar continuum emission. The result is a spuriously high F(4364) 
leading to the spuriously high values for T e . Further out, the temperatures show the red and 
yellow colour code. The pattern that looks like an “H” -shape region, or for the imaginative, 
the “batman symbol”, appeal's to be point symmetric with respect to the CS. This region 
is likely to be at the elevated T e s indicated because of the agreement in this central region 
with the STIS results shown in Figure 1. The full STIS long slit is shown superimposed 
on Figure 2. The tiles along the slit are entirely within the inner higher S/N region. The 
distances from the CS to the centres of the two symmetrical blue “lakes” that the STIS slit 
passes through are roughly 6.6". That distance corresponds to ~13 tiles on Figure 1, where 
the CS location is within tile 2G. Thus the troughs in the T, structure in Figure 1 near tile 
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13 and 39 correspond with the blue lakes. The temperatures also start to rise again further 
out as shown by the green colour. However, there is not a prominent increase in T e with 
increasing distance from the CS to the levels indicated from the T e distribution along the 
slit, in Figure 1. The sharp edges along the top and bottom of the inner elliptical bubble are 
blue and yellow in contrast to their neighboring mostly green pixels. This is likely due to a 
small image misregistration between the F437N and F502N images. The offset can occur for 
any number of reasons, including slightly different tilts of the filters in theii holdeis. The 
values for T e along the affected edges are too low at the top and too high at the bottom. 
The effect will be to raise the derived t\ slightly where these pixels are included. 

The determination of T e from I(5008)/I (4364) is not sensitive to the extinction correction 
for NGC 7009 because extinction is small. We may demonstrate this with an example using 
a r(H/3) of 0.1, which is close to the average value. Let us also use a case that results in 
a representative temperature 10,000 K. The observed F(5008)/F(4364) would be 153.543 to 
result in I (5008)/! (4364) = 148.057, which results in T e = 10,000 K by equation (8). If r(H0) 
were respectively 0.0 or 0.2 instead with the same observed flux ratio, the resulting T e would 
be 9891 and 10,1 12 K. Thus were we to carry out this analysis with a constant r( H/i) = 0.1, 
or no extinction correction at all, the results will be hardly altered. 


6. FRACTIONAL MEAN-SQUARE TEMPERATURE VARIATIONS 

Our STIS and WFPC2 analysis above present results in the plane of the sky. The 
observations hen; do not address temperature fluctuation along the line of sight, which may 
be characterized in terms of the average temperature To and fractional mean-square 7 ,■ 
variation ( t 2 ) as defined by Peimbert (1967). 


Tn = 


j T e N e N t dV 

f N'NidV 


,, _ r (T e -n)*N e Kdv (10) 

Tq f iV e A \dV 

where N, is the ion density N{() ++ ). The integration in equations (9) and (10) is over the 
column defined by each pixel (or tile for the STIS analysis), and along the line-of-sight (lot). 
We are unable to measure the i 1 along the los for any column (cross section 1 pixel or tile). 
If there are t 2 along the los, we can say that T(4364/5008) > 7), (e.g., Peimbert 1967; Rubin 
el, al. 1998). 

for each pixel/tile, we have calculated T(4364/5008). Then the intrinsic flux 1(5008), 
fully correcting F(5008) for extinction (see equ. 7), in each pixel/tile is used in conjunction 
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with 7; = T(4364/5008) for that pixel/tile, and assumed constant along the las, to derive 
the following: 

/(5008) = K (5008) j N e N t T~ 0 ' 5 exp(~x/ k T e ) dl = 7^(5008, T e ) j A',, N t dl . (11) 

Here x is the excitation energy above the ground state for the upper level of the 5008 
transition, k is Boltzmann’s constant, 7t (5008) is known from atomic data, and K (5008, 7 e ) 
has finally incorporated the known T e factor with the atomic constants. Here we again make 
the safe assumption of the low-N e limit (negligible collisional deexcitation) discussed eailiei. 

Let us now define analogous terms to represent the average 7 e ( I u,i ) and fractional 
mean-square T ( , variation (t\) in the plant; of the sky. 

_ f f T, N e N l dl dA ( 

M ff N e Ni dl dA ’ 

_ f f (Z - T 0:A f N e K dl dA f J TJ N, cH dA ] 

4 T* A ff N e Ni dl dA T* A ff N'NtdldA 

where dA represents an element of surface area in the plane of the sky and the integration 
over dl is for each pixel/tile along the los. 


6.1. WFPC2 Analysis 

For this part, the image in Figure 2 was further rotated to have the NGC 7009 major axis 
align with the x-axis. We obtained from this an array of intrinsic fluxes 1(4364) and 1(5008) 
for input to our own analysis code. For each of the above set of pixel-by-pixel solutions, we 
have the [O III] T e (that assumed no T e variations along the los). A circular region (radius 
of 10 pixels = 1 ") surrounding the CS is excluded from further calculations. Additionally, 
we exclude pixels near the periphery, where the S/N for the 4364 A line is poor. 

Calculations for To , a and t\ were performed using various distarrces front the CS ( ; . 
This was done several ways: assuming concentric circular shapes and concentric elliptical 
shapes, which is perhaps more realistic for this source. The bottom-line results are virtually 
unchanged between these different methods. In Figures 3 h 4, we present our results for 
To'A and (\ for an elliptical shape with an axial ratio of minor/ major = 0.808, matching 
our T,, map (Figure 2). The abscissa R. 0 U f. er = \/ a outer h ou i cr , where ci outer and b (iuter ate the 
semi-major and semi-minor axes. Out to the most distant elliptical contour used in our 
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analysis {a outeT = 14.5", 6 0U ,. er = 11.7", R olt(er = 13.0"), there are 52,995 good pixels with all 
(() III] T e values between 5531 16,382 K. 


The dashed lines pertain to values for T 0tA and t\ for the entire area interior to the 
ellipse with a uuter , b mUeT and is plotted at the single distance R wu<er (again excluding the CS) 
and represent cumulative results. The solid lines pertain to 70,. 4 and t A in the annular area 
between two adjacent ellipses separated by 5 pixels (0. 5 arcsec) along the major- 
axis. This is plotted as a single point at a distance R< X ntrmd — ( a outer b outeT + dinner b mne r 


) • 


There is nearly a rnonotonic decline in T 0<A with increasing distance (as characterized 
for the ellipses) from the CS. The results indicate very small values - < 0.012 foi t A 
throughout. The overall T M = 9912 K, i\ = 0.00360 for the 52,995 good pixels within the 
most distant elliptical contour used (excluding CS region). There is a steep climb in t' A in 
the outer elliptical annuli (solid line in Fig. 4), although this may be due to noisier data at 
larger distances from the CS. 


6.2. STIS Analysis 

Our STIS data allow an even more direct determination of T e and hence of 7o,a> 
albeit for a much smaller area than with VVFPC2. The average [O III] tempeiatuie To, a 
using 45 tiles (excluding the CS and fiducial bar) is 10,139 K; t 2 A is 0.00350. The median 
value (without regard to proper weighting) is 10,295 K. All 4 emission lines peak in tile 37, 
with tile 38 a close runner-up. This is where our slit crosses the “ridge” near the blue- 
green (He II 4687 [O ITT] 5008) boundary in the 3-colour image (Balick et, al. 1998) at 
roughly 5-6" from the CS. These tiles have relatively low [O III] T e of 9823 and 9755 K 
respectively (see Figure 1). No doubt, the enhancement in line emission here is due to a 
greater column density (and/or higher density) and not temperature. As determined from 
the extinction-corrected 5008 surface brightness (equation 11), these tiles have the highest 
(A r ( . dl ) -weighting in equations (12) and (13). 

We note that the single pixel (out of 52,995) with the highest 5008 surface brightness 
also has the highest (N t dl) -weighting, (T e N e N t dl) weighting, and (7/ N e N t dl)- 
weighting in the WFPC2 analysis. This pixel is located with respect to the CS 5.9" along 
the west major axis and 4.7" orthogonally below and is associated with the “ridge” near the 
blue green boundary of the Balick et al. (1998) picture. 
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7. DISCUSSION AND CONCLUSIONS 

We determined the electron temperature from the flux ratio of [O ITI] (4364/5008) using 
both STIS and VVFPC2 observations of NGC 7009. The resulting T e values from both the 
YVFPC2 map and the distribution along the STIS slit do not vary much, with almost all 
values between 9000 11,000 K. The higher T e values are closer to the central star and 

appear from WFPC2 images to coincide with the inner He ++ -zone. The temperature is 
higher in the He ++ -zone because heating is stronger due to radiation beyond 54.4 eV, while 
the cooling is less efficient as 0 ++ gives way to 0 +3 . Models show that some 0 ++ remains 
inside the He ++ zone, so the higher temperature there will be reflected in T(4364/5008) 
see, for example, figs. 11 and 14 of Harrington et al. (1982). 

The observations here do not address T e fluctuation along the line of sight. YVe assume 
for each square column (projection of 1 pixel for the WFPC2 data and 1 tile for the STIS 
data on the plane-of-the-sky) that the plasma along the line of sight is isothermal at the 
T(4304/5008). The analysis of both data sets for the average T e and fractional mean-square 
T e variations in the plane-of-the-sky, which we call T 0j a and t 2 A , gave consistent results. For 
the WFPC2 analysis, the preferred method used concentric elliptical contours to conform 
to the natural elliptical shape of the nebula. For Tq : a, the STIS average was 10,139 K, 
while the WFPC2 overall average to the outermost elliptical contour considered was 9912 K. 
Interior to the contours closer to the central star, l\j t A is higher as depicted in F igure 3. As 
described, in all the statistical analyses, a region around the CS is excluded. For both STIS 
and WFPC2, t\ is found to be very small, < 0.012 everywhere. Using STIS, we obtained 
0.0035, which is close to the WFPC2 overall value of 0.0036 within the most distant elliptical 
contour used. 

Because the continuum correction for the F437N data is a major source of uncertainty in 
the process of deriving T e from the WFPC2 observations, we attempted to assess this using 
another approach, Instead of estimating the continuum correction theoretically for F437N, 
we assume that STIS and WFPC2 are perfectly calibrated. By this we mean that the Tot al 
STIS emission/F437N in Table 1 should be unity (instead of 1.089). We then assume that 
the continuum emission contribution to the common WFPC2/STIS area SUM (see §3 and 
Table 1) be the remaining amount after accounting for the STIS 4364 and Hy emission 
as they would be observed passing through the F437N/telescope system. This reduces the 
Contm./4364 value for SUM in Table 1 to 0.1525 and alters the proportions that would 
comprise the F437N emission to be 0.77 : 0.12 : 0.11 for the 4364 line : continuum : Hy 
leakage. This treatment results in a modification of equation (4) with the factor 1.089 
becoming 1.0 and the factor 0.03218 becoming 0.02316. 

We carried out the subsequent analysis with just this change to equation (4). The 
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results were that the highest T e values were lowered several hundred degrees while the lowest 
7 ( values were raised several hundred degrees. Consequently, the effect was to deciease t A 
values somewhat. Repeating the elliptical annuli analysis, we found the VVFPC2 overall 
value of 0.00303 (decreased from 0.00360) within the most distant elliptical contour used. 
For the individual annuli, t\ < 0.01 everywhere. For T 0i/ 1, the WFPC‘2 overall average 
to the outermost elliptical contour considered increased slightly to 9940 K (from 9912 K). 
Kssentally, there were only minor changes to the curves depicted in Figures 3 and 4. We aie 
lead from the above computations to conclude that our results are robust and not sensitive 
in the situation considered here to the correction for continuum emission in the F437N data. 

Fluctuations in T e along the los are inevitable. We can make some comments about 
how our results might be adjusted by T e variations along the los. The relationship between 
7" (4304/5008) and 7’ 0 is 


T(4364/5008) = T 0 [1 + 0.5(91, 200/T o - 3)f' 2 ] , (14) 

(e.g.. Peimbert 1907; Rubin 1969). For a given column, let us take a representative 
7'(4364 /5008) of 10,100 K, which is close to the average T U<A for the ST1S analysis. If 
7ji were 9000 K, t hen t 2 = 0.0343 along that los in the column; if T 0 were 8000 K, then t 2 = 
0.0625. If we assume that generally a given los column will be characterized by 7y = 9000 K, 
1100 K lower than 7^4364/5008) , then according to equation (13), one would expect t\ to 
| H . a factor of 1 .26 larger. This is due to a simple temperature scaling by a factor of T 2 A /T$ . 
Our conclusion will still remain that is very small as measured by the FIST data here. 


We do not have the data here to characterize T e variations in 3-dimensions. One way 
that this may be realized is to consider the following case. If every single pixel/tile los- 
column were to have T(4364/5008) = 10,100 K, T 0 = 9000 K, and t 2 = 0.0343, the analysis 
for variations in the plane-of-the-skv that we present here will result in T£ a = 10,100 K and 
t' 2 A : o. However, the overall true 3-dimensional T 0 = 9000 K and t 2 = 0.0343. 


Irrespective of the above contrived case, we draw the following conclusion. Although 
we have measured l l A in 2-dimensions and not yet t 2 in 3-dimensions, our results indicate 
that temperature fluctuations alone are unlikely to explain for NGC 7009 the fat ten of 
discrepancy (Liu et al. 1995) between heavy element abundances inferred from emission lines 
that are eollisionally excited compared with those that are due to recombination/cascading. 
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uescnption ol EIXHUNTER 

This was a 2001 summer project of Aaron Svoboda working at NASA/Ames with RR. 
lere are separate procedures to identify and remove these pixels from emission lines and 
rom the continuum For the continuum regions of interest, generally adjacent to emission 
hnes t, 0 be measured, sections of the rectified 2-D image are identified that are reasonably 

• , ” UX l ' lhe Spatml and dis l ,ersio11 directions. Sections are necessary in order to 
avoid the CS and ST1S fiducial bars. PIXHUNTER works by computing the mean pixel 

t ? T' U deV , iati0n ° fagiVen SeCtl ° n - After COII iputing these values, it again goes 

hiough the image and outputs to a file the coordinates of any pixel whose flux value is 

fuither away from the mean than ± a specified multiple of the standard deviation (a); here 
we used o a. We linearly interpolate the good data to replace the bad pixel values. 

Though the above algorithm for finding the bad pixels works well for the continuum 
see ions it . is incompatible with emission lines, because the flux in an emission line can vary 

PIXH I J NTFR 111 f WaVeIe, ' gth aild Spatia ' dimeustons - ^ cleaning an emission line. 
HUNTER works interactively with splat We plot and compare a single column (fixed 

wavelength) in the line at a time with one of the fitting functions available with splat (e g a 

Hm d!l?i. P< ; " r n F ° r r h in<JividUal C ° ,Um “ in «». this fitting function represents 

' . 0,1 ° f f UX W,th Spat,al l ,ositio11 alo »g the slit. Following the fit to the emission 

line distribution, splat will flag all pixels that deviate by more than a specified multiple of a 

'om the functional fit. Various keys/commands within splat permit further processing such 
as inter pointing all of the flagged pixels in a column to the function. When the resulting 
■ mage is satis acton ly repaired, it can lie saved back into the original 2-D rectified image 
Kach column can be treated in this fashion by scrolling through interactively using existing 
op aons/toos within splat. There is another method for cleaning bad pixels in the vicinity 
of pronounced peaks m the spatial distribution. This involves taking the ratio of a single 

C ° 1U, 7 t ,C hr * t: U> blkav U of the entire line. Again, splat is used interactively, operating 
on each such ratio, m a similar fashion to the first method. 
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